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The Ashtekar-Krishnan energy-balance law for dynamical horizons, expressing the increase in 
mass-energy of a general black hole in terms of the infalling matter and gravitational radiation, is 
expressed in terms of general trapping horizons, allowing the inclusion of null (isolated) horizons as 
well as spatial (dynamical) horizons. This first law of black-hole dynamics is given in differential 
and integral forms, regular in the null limit. An effective gravitational-radiation energy tensor is 
obtained, providing measures of both ingoing and outgoing, transverse and longitudinal gravitational 
. . . radiation on and near a black hole. Corresponding energy-tensor forms of the first law involve a 

, preferred time vector which plays the role for dynamical black holes which the stationary Killing 

' vector plays for stationary black holes. Identifying an energy flux, vanishing if and only if the horizon 

' is null, allows a division into energy-supply and work terms, as in the first law of thermodynamics. 

The energy supply can be expressed in terms of area increase and a newly defined surface gravity, 
^ ■ yielding a Gibbs-like equation, with a similar form to the so-called first law for stationary black 

Q ' holes. A Clausius-like relation suggests a definition of geometric entropy flux. Taking entropy as 

\ area/4 for dynamical black holes, it is shown that geometric entropy is conserved: the entropy 

. of the black hole equals the geometric entropy supplied by the infalling matter and gravitational 

^SJ ' radiation. The area or entropy of a dynamical horizon increases by the so-called second law, not 

OA ' because entropy is produced, but because black holes classically are perfect absorbers. 

CN . PACS numbers: 04.70.Bw, 04.30.Db, 04.70.Dy 
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' Black holes are perhaps the most exotic and energetic objects in the universe. Their theoretical history is long and 
. winding: the earliest such solution to the field equations of Einstein's theory of General Relativity 1] was found by 
Schwarzschild 2] almost immediately, but not understood as such for decades 3]. In a few years around 1970, there 
Q was rapid theoretical progress, with the introduction of the term black hole by Wheeler [J| and the development of 
the classical four laws of black-hole mechanics [HSSBj supposedly analogous to the laws of thermodynamics. Since 
^' ^ ^ then, astrophysical evidence has increasingly accumulated not only for stellar-mass supernova-remnant black holes, 
bi)' but for supermassive black holes, mysteriously present at the heart of most if not all galaxies and powering active 
■ galactic nuclei . Cataclysmic events such as binary black-hole mergers are predicted to produce gravitational waves 
• ^ , which are observable on or near our home planet, for which a new generation of detectors is being developed |lCll |. 

Consequently, recent years have seen a great deal of work on numerical simulations to study how black holes evolve 
^ , according to given initial conditions, and what gravitational radiation they may produce (llj . 
■ - - ■ Such progress leaves the textbook theory of black holes seriously out of date. Much is known about stationary black 
holes, for instance the zeroth and first laws just mentioned, but dynamical black holes are much more complex. Of 
the classical laws, only Hawking's area theorem has generality, but it applies to event horizons, which are theoretical 
constructs which cannot be located by mortals. It is quite timely that Hawking has recently recanted, writing that 
"a true event horizon never forms, just an apparent horizon" |l2l |. Unfortunately, Hawking's definition of apparent 
horizon is also not the most appropriate to define black holes, due to its global nature and slicing dependence; for 
instance, the Schwarzschild black hole may be globally sliced so that there is no apparent horizon [l^]. 

About ten years ago, the author began a program to understand local, dynamical properties of black holes [T3.IT5l| . 
The basic idea is that black holes contain trapped surfaces, where both ingoing and outgoing light wavefronts are 
converging, and that one can locate the surface of the black hole by marginal surfaces, where outgoing light rays 
are instantaneously parallel. A trapping horizon is a hypersurface foliated by marginal surfaces. Locally classifying 
trapping horizons as future or past, and outer or inner, it was proposed that a future outer trapping horizon char- 
acterizes non-degenerate black holes. Some general results were that: there are future trapped surfaces just inside 
such a horizon; the horizon is achronal, being null only in the locally stationary case and otherwise spatial, assuming 
the null energy condition; the marginal surfaces have spherical topology, assuming the dominant energy condition; 
and the area A of the horizon is non-decreasing, A' > 0, and increasing if spatial. The last property is analogous to 
Hawking's area theorem, but for a practically locatable horizon. Trapping horizons can be numerically located by 
so-called apparent-horizon finders |l6l Il7|. which actually find marginal surfaces; they do not check every surface in 
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the hypersurface to see whether it is outer trapped, as required by the definition of apparent horizon. Under smooth- 
ness assumptions 0,^3, apparent horizons are marginal surfaces, but not vice versa. Incidentally, the resolution to 
the supposed black-hole-information paradox is simple, using versions of the above results for matter violating the 
null energy condition: as a black hole evaporates, the ingoing negative-energy Hawking radiation causes the trapping 
horizon to shrink and become temporal, so that information can cross it in both directions |il9„ .2(1] . There never was 
a paradox, just a fundamental misunderstanding, that black holes are usefully defined by event horizons. 

In terms of trapping horizons, a comprehensive picture of black-hole dynamics was first developed in spherical 
syinmetry. In this case, there are local definitions of active gravitational mass-energy E (21.] and surface gravity k 
|2^ which have many physically expected properties. The gradient dE of energy, expressed in terms of the energy- 
momentum of the matter, divides naturally into an energy-supply term Atp and a work term wdV {V = ^ttR^, 
A = AttR'^) such that the energy fiux tp vanishes on a trapping horizon if and only if it is null. This energy-balance 
equation was called the unified first law for various reasons: projecting it along the flow of a thermodynamic fluid 
yields a first law of relativistic thermodynamics; projecting it along null infinity yields the Bondi energy-loss equation, 
with ip reducing to the Bondi flux; and projecting tp along a trapping horizon gives kA' /Sir, yielding an equation 
E' — kA' /8tt + wV with the same form as the so-called first law for stationary black holes, which is really analogous 
to the Gibbs equation rather than the first law of thermodynamics. Including a zeroth law the possible local 
properties of dynamical black holes, independent of particular matter models, are thereby known. This can also 
be achieved in cylindrical symmetry |24j . where the energy flux ^ divides into contributions from the matter and 
the gravitational radiation. A quasi-spherical approximation 25^ _2fii 21i .28] also allows generalizations of all these 
geometrical and physical properties of black holes and gravitational radiation. An effective energy tensor for the 
gravitational radiation can be given in all these cases. 

Ashtekar & Krishnan recently found an energy-balance law for dynamical horizons p9ll30ll3lll32j ]. defined as spatial 
future trapping horizons. This arose from earlier work on isolated horizons, types of null trapping horizon for which 
generalizations of the classical laws of black-hole statics were found [33, IM, H^, Hg, H^, HI] . The new energy-balance 
law describes how a general black hole grows due to the infalling matter and gravitational radiation. In this article, 
these two threads are drawn together, in particular deriving the first law for completely general trapping horizons, 
so as to include both spatial (dynamical) and null (isolated) horizons, as well as horizons of white holes, traversable 
wormholes, cosmological models and evaporating black holes (§VI). No displayed equation will assume any restriction 
on the type of trapping horizon. 

A terminological mismatch should be mentioned at the outset: the first law here is what Ashtekar & Krishnan 
called a balance equation for area or energy, while their first law generalizes the so-called first law of black-hole 
mechanics, involving new definitions of angular momentum and surface gravity, which are not considered here. They 
gave an integral form using proper volume, which becomes singular in the null limit, so here the first law is written 
in differential form and in an alternative integral form, which are both regular in the null limit (§VII). An effective 
gravitational-radiation energy tensor Q and a preferred time vector x are obtained, yielding energy-tensor forms of 
the first law (§VIII). A division into energy-supply and work terms, generalizing the above structure in spherical or 
cylindrical symmetry, allows a Gibbs- like equation involving a new definition of surface gravity k (§IX). The energy 
fiux of the matter and gravitational radiation satisfies a Clausius-like relation involving /t, suggesting a definition 
of geometric entropy flirx 2-K'tp/ k. Then it is found that geometric entropy is conserved: the geometric entropy of the 
black hole equals the geometric entropy supplied to the black hole by the infalling matter and gravitational radiation 
(§X). The results are preceded by brief reviews of basic thermodynamics (§11), the employed formalism of dual- null 
dynamics (§111), the definition of trapping horizon (§IV) and the area and signature laws (§V), and followed by a 
Conclusion (§XI). See j39| for a short version and 113,143] for different but related approaches. Standard Einstein 
gravity is assumed, though the ideas generalize. 



II. BASIC THERMODYNAMICS 



It seems appropriate to begin with a brief summary of basic thermodynamics, due to the parallels often drawn for 
black holes, and the fact that they are not always accurate. See fuller treatments [43. l43l 143. lisi] and beware any 
source which formulates the laws of thermodynamics using state-space differentials d and meaningless derivatives 5. 

In classical thermodynamics, the basic quantities are temperature heat supply Q, work internal energy 
(actually thermal energy or simply heat) H and entropy S. The classical first law is 

H = Q + W (1) 

where the dot denotes the material or comoving derivative. For instance, for an inviscid fluid, the work is given by 
W = —pV, where V is the volume and p the pressure of the fluid, so that the first law reads 



H = Q-pV. 



(2) 
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The classical second law, originally due to Clausius, who used it to define entropy, is 

S > Q/i?. (3) 

These integral forms of the laws respectively require the pressure and temperature to be spatially constant. 
The entropy may be divided into entropy supply So, given by 

- O/i? (4) 

and entropy production S ^ So- Then the second law may be written as 

^ > 5o (5) 

which expresses entropy production. In words: S is the entropy of the system, where system means a comoving volume 
of material, and So is the entropy supplied to the system. Thus the second law implies that the total entropy of an 
isolated system, such as the whole universe, cannot decrease. Here it should be stressed that dynamical black holes are 
not isolated systems, since they absorb energy and entropy. Then the property that black holes have non-decreasing 
area, A' > 0, normally called the second law of black-hole mechanics, is actually not analogous to the second law of 
thermodynamics. Entropy production and entropy increase have entirely different meanings for non-isolated systems. 

Equality in the second law holds in thermostatics, traditionally called equilibrium thermodynamics or reversible 
thermodynamics. In the thermostatic case, the first and second laws for an inviscid fluid imply 

H = dS-pV (6) 

which is the Gibbs equation, or rather its material or comoving form. Note that what is normally called the first 
law of black- hole mechanics for stationary black holes 5], involving area A — AS and surface gravity k, = 27ri9, is 
actually analogous to the Gibbs equation, rather than the first law of thermodynamics. The latter does not involve 
temperature or entropy, but simply expresses energy balance. 

Thermodynamics can be formulated as a local field theory, with H and S replaced by thermal energy density and 
entropy density respectively, Q replaced by a thermal flux vector q such that 



Q = - j) *n- q (7) 
and So replaced by an entropy fiux vector ip — q/-d such that 

So = — f *n ■ ip (8) 



where the integrals are over a surface bounding the system, with vector area element *n. In terms of these and other 
fields, e.g. density, velocity and stress for a fluid, the first and second laws and the Gibbs equation can be localized 
pL^s .Ml 1^ ' These three localized equations or inequalities can then be used to derive dissipative relations, in 
the simplest case the Fourier equation for q and the Newtonian-fluid equation for the viscous stress, leading to the 
Navier-Stokes equation. Thus it should be stressed that the first law and comoving Gibbs equation are still assumed 
fundamentally and fruitfully in true (non-equilibrium) thermodynamics as well as in thermostatics. Widespread 
folklore to the contrary is sometimes used to argue that the so-called first law of black-hole mechanics, obtained as a 
property of stationary black holes, should not be expected to generalize to dynamical black holes. Again, this does 
not constitute a correct analogy with true thermodynamics. 



III. DUAL-NULL DYNAMICS 



Trapping horizons are generally defined as hypersurfaces which may have any causal nature, foliated by marginal 
surfaces. To study them, it is useful to employ the formalism of dual- null dynamics 0jli3j describing two families 
of null hypersurfaces, intersecting in a two-parameter family of transverse spatial surfaces, as summarized in this 
section. There are various reasons: marginal surfaces are defined as extremal surfaces of null hypersurfaces; a spatial 
trapping horizon locally determines a unique dual-null foliation, generated from the marginal surfaces in the null 
normal directions; and the null limit, where a dynamical horizon reduces to an isolated horizon, is naturally included 
in the formalism, whereas more conventional treatments of spatial hypersurfaces become degenerate in the null limit, 
basically because normal vectors become tangent. For a null trapping horizon, the dual-null foliation is not unique, 
so subtleties remain in describing partial spatial, partially null trapping horizons. 
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Denoting the space-time metric by g and labelling the null hypersurfaces by coordinates which increase to the 
future, the normal 1-forms 



n 



dx^ (9) 



therefore satisfy 

g-\n^,n^) = Q. (10) 

The relative normalization of the null normals may be encoded in a function / defined by 

ef = -g-\n+ ,n-) (11) 

where the metric sign convention is that spatial metrics are positive definite. Some readers may prefer to write 
(7"' = g~^(n+,n~) for more manifest invariance and remember that g'^ < 0. The induced metric on the transverse 
surfaces, the spatial surfaces of intersection, is found to be 

h = g + 2e~fn'^ ® n~ (12) 

where ® denotes the symmetric tensor product. The dynamics are generated by two commuting evolution vectors 
u±: 

[w+,M_]=0 (13) 

where the brackets denote the Lie bracket or commutator. Thus there is an integrable evolution space spanned by 
(u-i-, u_). There are two shift vectors 

s± = _Lu± (14) 
where _L indicates projection by h. The null normal vectors 

l± = u± — s± = g~^{n^) (15) 

are future-null and satisfy 

9{l±,l±) = (16) 

g{l+,l.) = -e-f (17) 

l±-dx^ = 1 (18) 

l±-dX^ = (19) 

_LZ± = (20) 

where a dot denotes symmetric contraction. In a coordinate basis (u+,u_,ea) such that u± = d/dx"^, where Ba = 
d/dx"" is a basis for the transverse surfaces, the metric takes the form 

g = habidx" + s'^dx+ + s"Ldx~) (g) {dx^ + s''_^_dx+ + s^_dx~) - 2e~^ dx+ O dx~. (21) 

Then (h, /, s±) arc configuration fields and the independent momentum fields are found to be linear combinations of 
the following transverse tensors: 

6l± = *L±*1 (22) 
a± = ±L±h-9±h (23) 
u± = L±f (24) 
u = ie//i([/_,Z+]) (25) 

where * is the Hodge operator of h and L± is shorthand for the Lie derivative along l± . Then the functions 6± are the 

expansions, the traceless bilinear forms a± are the shears, the 1-form lo is the twist, measuring the lack of integrability 
of the normal space, and the functions v± are the inaffinities, measuring the failure of the null normals to be affine. 
The fields {0±, a±, i^±, w) encode the extrinsic curvature of the dual- null foliation. These extrinsic fields are unique up 
to duality ± i-^ =f and diffeomorphisms x^(x^) which relabel the null hypersurfaces. It will also be convenient 

to use the normal fundamental forms 

C(±)=e-^±((nT-V«)n±)=D//2Tc^ (26) 
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FIG. 1: A dual-null foliation: the commuting evolution vectors u± — d/dx"^ generate the transverse surfaces 5", while their 
null normal projections l± generally do not commute. 



where V and D are the covariant derivatives of g and h respectively, and V A rv^ ~ has been used. Here (±) 
indicates a label, not an index; Qj-) are two generally distinct transverse 1-forms. One can compose them as a 2-form 
(3 in the normal space, defined by (3{jjL,v) = ■ V)i^) for normal 1-forms {^jv), with components (3^^ = c-'^Cct)' 

= 0, but such notation becomes cumbersome. Likewise, one can compose the expansions and shears into a second 
fundamental form, but it is more convenient to separate them. 

One subtle point concerns the evolution vectors u± versus the the null normal vectors l±, differing by the shift 
vectors s±. In a numerical evolution, one would be evolving using the field equations with -LLu^ derivatives on the 
left-hand side, since such Lie propagation of a point takes it to other points with the same angular coordinates. In 
particular, evolving along then u_ takes one to the same point as evolving along m_ then m+, since they commute, 
as depicted in Fig.^ However, l± generally do not commute (as measured by uj), so that evolving along Z+ then l- 
takes one to a generally different point, though in the same transverse surface, as evolving along Z_ then l+. On the 
other hand, for analytical purposes it is easier to use l± than u±, writing the field equations with -LL± derivatives 
on the left-hand side. The same issue exists in the 3-f-l formalism, where the evolution vector differs by the lapse 
function and shift vector from the unit normal vector. Another subtle point is that l± are not general tetrad vectors, 
since they are defined in terms of n^, which must be closed. In particular, this means that / cannot be fixed to zero 
for a general dual-null foliation. Another way to see this is that its derivatives L±L^f are determined by the Einstein 
system in terms of the free initial data , even in spherical symmetry. 

The dual-null Hamilton equations and integrability conditions for vacuum Einstein gravity were derived previously 
[47| |. with matter terms added subsequently |27l| . Denoting projections of the energy tensor T by T±± — T(l±,l±) 
and = T(l^, /_), the relevant components of the field equations are just 

L±9± = -i^±0±-0l/2-\\a±\\yA~87rT±± (27) 
L^0± = -0+0- - e-f {n/2 - |C(±) I' + D ■ Cf±)) + 87rT+_ (28) 

where is the Ricci scalar of h (conventionally positive for spheres), a sharp (jj) denotes the contravariant dual with 
respect to = (index raising), |CP = C • ^-nd ||(t|P = cr : cr", where the colon denotes double symmetric 
contraction. Units are such that Newton's gravitational constant is unity. The first equation is the well-known null 
focusing equation and the second has been called the cross-focusing equation 0, 0| . The null energy condition 
implies 

T±± > (29) 
and the dominant energy condition additionally implies 



r+_ > 0. 



(30) 
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IV. TRAPPING HORIZONS 



The dual-null formalism may be applied to one-parameter families of transverse surfaces generated by a vector 

^^ei++ri- D^^^O. (31) 

This means that ^ = d/dx is normal to the constant-a; transverse surfaces, so that can be taken outside transverse 
surface integrals ^. The area of the transverse surfaces is 



(h *1 (32) 



and the area radius 



is often more convenient. The Hawking energy [4^ 



R ^ ^/A/Att (33) 



E I ^(^a + e^e+e) (34) 

lOTT J 

will be used as a measure of the active gravitational mass-energy on a transverse surface. On a stationary black-hole 
horizon, it is also known as the irreducible mass: the mass which must remain even if rotational or electrical energy 
is extracted. 

Consider a trapping horizon generated by a vector — d/dx, so that the constant- a; surfaces are marginal surfaces, 
where one of the null expansions 9± vanishes. This leaves the freedom to relabel the marginal surfaces, x i— > x{x), 
under which all the key equations will be manifestly invariant. Equations holding on a trapping horizon will be 
denoted by the weak equality symbol =. Initially, the case 6*+ = will be considered in detail, with the case 6*- = 
included subsequently. The fundamental equation describing the evolution of a trapping horizon is 

{)'^L^e+ = i+L+e+ + cL-e+ (35) 

where denotes the Lie derivative along ^. This will be used together with the Einstein equation to derive the first 
law for any trapping horizon. Even without the Einstein equation, one can use it to deduce relationships between the 
signs of and L±9^, which determine the causal nature of the trapping horizon and whether its area increases or 
decreases, via 

L^A^ j *{te++CO-). (36) 

For clarity, all such inequalities are collected in the next section, so that the remainder of the article applies to any 
trapping horizon. 



V. AREA AND SIGNATURE LAWS 

Trapping horizons were previously classified 'l4l into one of four non-degenerate types: future (respectively past) 
if 6'- < (respectively 0- > 0) and outer (respectively inner) if L_0+ < (respectively L^9+ > 0) on the trapping 
horizon = 0. For each type, there are trapped surfaces {0+9- > 0) to one side of the horizon and untrapped surfaces 
{9+9- < 0) to the other side, which is not guaranteed if the inequalities are relaxed even to non-strict inequalities. 
The causal type of the horizon is determined pointwise by the relative signs of : spatial if they have opposite signs, 
null if one vanishes and the other does not, and temporal if they have the same (non-zero) sign. Since the null energy 
condition H29|l and focusing equation (|27|l imply L+9+ < on the trapping horizon, it follows from the fundamental 
equation (|35|l that outer trapping horizons are achronal (spatial or null), while inner trapping horizons are causal 
(temporal or null); the signature law Furthermore, fixing the orientation of C by C"*" > Oj it follows from l|3t)|) that 
the area of a future outer or past inner trapping horizon is non-decreasing, L^A > 0, while the area of a past outer or 
future inner trapping horizon is non-increasing, L^A < 0; the area law As corollaries, the horizon is null and has 
instantaneously constant area if and only if the ingoing energy density + vanishes, where &++ = ||(T+|p/327r 
(|68|l can be understood subsequently as the effective energy density of ingoing gravitational radiation. 

As mentioned in the Introduction, non-degenerate black holes may be characterized by future outer trapping 
horizons. Ashtekar & Krishnan instead defined dynamical horizons as spatial future trapping horizons. Then 0_ < 
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on the horizon and ^ have opposite signs. Choosing the orientation of l± such that Z+ is outward and Z_ inward, 
^+ > 0, ^~ < and it follows directly from (|36l) that the area is increasing, L^A > 0, a strict version of the above 
area law. Actually, for black holes, one is normally interested in future trapping horizons which are either spatial 
(dynamical) or future- null (isolated), or partially spatial and partially null. In such cases > and < 0, which 
immediately gives the non-strict area law L^A > 0. Note also from the fundamental equation (|35|l that a dynamical 
horizon satisfies < under the null energy condition, so it is either a future outer trapping horizon or degenerate. 

The degenerate cases allow dynamical horizons in space-times without trapped surfaces 49], reflecting the need for 
something like the outer condition to characterize a black hole. In practice, the outer horizon of a black hole is likely 
to satisfy both definitions, except when it becomes stationary or instantaneously stationary. 

An evaporating black hole may also be described using trapping horizons. The only difference with the above 
discussion is that the null energy condition is violated by Hawking radiation, for which the ingoing radiation has 
negative energy density, T-|_+ < 0. Assuming that this dominates the positive energy density of ingoing gravitational 
radiation, T_|__|_ -|- &++ < 0, the focusing equation H27() implies > on the trapping horizon. For an outer 

horizon, the fundamental equation (|35|l implies that have the same sign, so that the horizon is temporal, while 
(|36|1 shows that the area is decreasing for a future horizon, L^A < 0. Thus the black-hole horizon is shrinking and 
two-way traversable. Clearly matter can escape from an evaporating black hole. The strange belief that information 
cannot escape from an evaporating black hole seems to be based on the impractical event-horizon definition of black 
hole as a region of no escape. 



VI. FIRST LAW: ENERGY FLUX AND WORK 



Henceforth completely general trapping horizons will be considered, so that all the following displayed equations will 
apply not only to outer black-hole horizons under the usual energy conditions, but to inner black-hole horizons, white 
holes, cosmological horizons, wormholc mouths and evaporating black holes. Expanding the fundamental relation 
using the focusing equations l(77H^ yields 

- L^9+ = -e+ (87rT++ + ||<t+||V4) + T (87rr+- - e'^ (5R/2 - |Cp + D ■ C^)) (37) 

where ( — temporarily simplifies the notation. Multiplying by c^/Stt and integrating over the transverse surfaces, 
using the Gauss-Bonnet theorem ^ *'St — Sir, the Gauss divergence theorem ^ *D ■ a = and rearranging yields 




Since = in the null case, this shows that both and (t+ must vanish on a null horizon 14], assuming the null 
energy condition. Here spherical topology has been assumed; otherwise, for compact oricntable transverse surfaces, 
the right-hand side of (|38|l is multiplied by 1 — 7, where 7 is the genus or number of handles. Then the dominant energy 
condition implies 7 < 1, leading to the topology law [lj|: the transverse surfaces are either spherical or toroidal, the 
latter case requiring very special conditions (including vanishing Gaussian curvature) and anyway being excluded for 
(non-degenerate) outer trapping horizons, < 0. 

The Hawking mass-energy (|34f) satisfies 

E = R/2 (39) 

on a trapping horizon, which can be regarded as a generalization of irreducible mass from stationary to non-stationary 
black holes, since the area law ensures its irreducibility under the null energy condition. Then the identity (|38|l yields 

This is a dual- null differential version of the energy-balance law found by Ashtekar & Krishnan |29|, |30j , as compared 
more explicitly in the next section. One may fix the normalization / = and, for a spatial horizon, one may fix the 
scaling of the null normals such that S,'^ = —1 and the generating vector ^ such that L^R = 1. In the following, 
all gauge freedom will be retained for generality, but readers may wish on a first reading to mentally set / = and, 
if interested only in spatial trapping horizons, ^~'~/^~ — — 1- The four terms are all geometrical invariants of the 
dual-null foliation, as shown explicitly below, and therefore of the horizon (as an embedded hypersurface) unless it 
becomes null. Since the formalism is manifestly covariant on the transverse surfaces, checking invariance reduces to 
writing e-'^ = — g"' and matching ± indices. 



(40) 
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The four terms in parentheses in H40() are all manifestly positive, assuming the dominant energy condition. The 
term gives the energy flux of the matter propagating in the null direction into the horizon. Consequently it is natural 
to interpret the (t+ shear term as giving the energy flux of the transverse gravitational radiation propagating in the 
null direction into the horizon. This term has the same form as that of the Bondi flux of gravitational radiation at null 
infinity 50., 51] , the same form as a localized energy flux of gravitational radiation in a quasi-spherical approximation 

f, and a similar form to the energy flux of linearized gravitational radiation in the high-frequency approximation 
so its physical interpretation seems sound. The term gives a matter energy density, so the C term can be 

interpreted as giving a corresponding gravitational energy density. Ashtekar & Krishnan interpreted it as also due 
to gravitational radiation, and here it will be interpreted as an energy density of longitudinal gravitational radiation. 
This is much less familiar than transverse gravitational radiation and is absent in all the above approximations, but 
the interpretation can be understood in a spin-coefficient formulation, to be presented elsewhere 52] . It should be 
mentioned that there is a widespread belief that longitudinal gravitational radiation does not exist in Einstein gravity, 
apparently due to an argument in linearized theory that the longitudinal modes are purely gauge-dependent, and 
the fact that only the transverse mode contributes to the Bondi flux. On the first point, Szekeres characterized 
ingoing and outgoing, transverse and longitudinal gravitational radiation by their effect via the geodesic equation on 
a "gravitational compass" of test particles jMj,Mj- On the second point, it can be shown that the energy densities of 
the outgoing transverse and longitudinal modes fall off near future null infinity as and l/R^ respectively |^ [s^] ■ 
The expression (|40ll separates the first term in parentheses, which vanishes for null horizons (assuming the null 
energy condition), from the second term in parentheses, which is generally non-zero for horizons of any causal nature. 
This separation need not appear for spatial trapping horizons, but it will be stressed in the following, since the null 
case is a physically important limit, where dynamical horizons reduce to isolated horizons, or more prosaically, where 
a growing black hole ceases to grow. 

The next task is to write the new law in a more manifestly invariant form. The spherically symmetric case is a 
useful guide; there the unified first law was found as dE = Aip + wAdR for certain invariants ■0 and w of the matter 
energy tensor . The corresponding formulae read 

ui„ = -trace T/2 (41) 
V^™ = T • {dRf + w„,dR (42) 

where the trace is in the normal space and the subscript m is introduced to refer to the matter. In the current 
generalized context, these quantities are invariants of the dual- null foliation. Explicitly, 

= e-^T+_ (43) 
(V'„)± = -efT±±L^R (44) 

where ip = ip^dx~^ + 'ip^dx~ . Comparing with the first law (|40|l . one can define corresponding quantities for the 
gravitational radiation by 



As in the spherically symmetric case [22j , one may call 

W Wm + Wg (47) 

the work density and 

i; = V^m + (48) 

the energy flux 1-form. 

On a trapping horizon with 0^ = 0, it follows that L+i? = 0, ip- = 0, ^ • i/' — ^nd L^R = S^~L^R, yielding 

^^^-|^e/fT,, + %^Ve^■ (49) 



327r 



Thus the first law 140|) becomes 



j) ■ t/j + <j) *wL^R 



L^E = 6 ■ iIj + 6 *wL^R. (50) 
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This is the first law of black-hole dynamics, in the desired geometrically invariant form. (The Lie derivative acting 
on integral scalars like E and R is just the partial derivative d/dx). The energy flux ^•■0 vanishes if the horizon is null, 
while the work density w is generally non-zero for horizons of any causal nature. The two terms § *^-ip and / *wL^R 
may be called respectively the energy-supply and work terms, by analogy with the first law of thermodynamics JQ). 

The above derivation of the first law applied to a trapping horizon with 0+ = 0. For a trapping horizon with 9- = 0, 
one obtains the same formula with a different w. For completeness one can define 

= e/T+_ + (51) 

and use w = W(^±) for horizons with 9± = 0. Recall that one can anyway fix the normalization / = 0, in which case 
uj = = — On the other hand, the same ^ appears in both cases, indicating that both components 'ip± are 
correct, i.e. that the energy densities given by the first bracket in (|40() . which at first sight are just scalars, can be 
naturally derived from an energy flux 1-form ip as ip- 



VII. INTEGRAL FORMS OF THE FIRST LAW 



Ashtekar & Krishnan derived an integral form of the first law, using proper volume, whereas the first law H5()(l 
derived above is in differential form. For comparison, it can be written simply as 



<j)*eL^R (52) 



where 

^ LfR 



w (53) 



is the combined energy density, where the division into energy-supply and work terms can be ignored in this section. It 
can be independently divided into matter and gravitational-radiation terms, e -t- , in the obvious way, yielding 
the explicit expressions 

em = ef (t+^ - ^T±±) (54) 



i: 

327r - 



^ ^|C(±)P-e^|^||a±|n (55) 



for a trapping horizon with 9± = 0. 

A corresponding integral form of the first law is 



[E] J *€L^R A dx (56) 



which expresses the change in E along the horizon, from one marginal surface to another. In this article, J always 
denotes such a hypersurface integral between transverse surfaces, f always denotes a transverse surface integral and 
[] denotes the change in such a quantity between transverse surfaces. This manifestly invariant expression uses the 
generator- volume element *1 A dx. Alternatively, if one wishes to use the proper- volume element 

*1 = *1 A ds = *y/gxx A dx (57) 

where s is arc length along the horizon- generating vector ^ — d/dx and 



V9^ = VgKO = ds/dx (58) 
is the corresponding scale factor, then the integral first law can be written as 

[E] = J a (59) 

where 

e = eLjjR = rj ■ ip + wLj^R (60) 
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is the proper energy density and 

V = (./Vg^ ^ d/ds (61) 

is the vector which is parallel to a normal generating vector and differentiates with respect to arc length, rj ■ ds — 1. 
Note that SI, s and r] are all independent of relabellings of the marginal surfaces, x ^ x{x). 

The last form is a compact version of the Ashtekar-Krishnan energy-balance law, so one can finally check consistency. 
They assumed gauge choices which here correspond to / = 0, C^/C^ = —1 and L^^R = 1, mentioned after H40|) 
as admissible for a spatial horizon, and their null normals are \/2^±, but it is now straightforward to check their 
expressions against the explicit expressions (|3Hj^ . which reduce to — and eg ^ (||(T+|p-|-4|Q_|_)p)/327r 

for = 0, with e reducing to e/y^gxx- They also gave a more general form for [a;/2], not fixing L^R, discussing 
a; = i? as a special case; this can easily be reproduced here simply by dividing both sides of the first law (in whatever 
form) by L^R. Note some notational pitfalls: their lapse N is \l ^gxx here, their coordinate r is a; here, their u is 
cr+ / \/2 here, and their permissible vectors ^ would here be Z+ ^2/ gxx, rescaling the outward null normal rather 
than the horizon-generating vector ^ (to ry = S,/ ^Jgxx)- Either way, the rescaled vector is ill-defined if the trapping 
horizon becomes null, gxx — > 0. In this limit, the proper-volume element vanishes, while the apparently ill-defined 
proper energy density e turns out to be finite. This suggests using either the generator-volume or differential form to 
deal with partially spatial, partially null trapping horizons. 



VIII. ENERGY-TENSOR FORMS: EFFECTIVE GRAVITATIONAL-RADIATION ENERGY TENSOR 

Dividing the integrated energy fiuxes into those due to the matter and the gravitational radiation in the obvious 
way, 

[E] = [EU + [E]g. (62) 

Ashtekar & Krishnan stressed that one could obtain the integrated matter flux as an equation which here would be 
[E]m = J *T{l^y^2/gxx, t) in their gauge choice, where f is the unit normal vector to a spatial trapping horizon. 
However, unit f does not exist for null trapping horizons. Nevertheless a natural normal does exist for any trapping 
horizon, namely the vector r dual to the generating vector ^ in the normal space: 

T = i^ (63) 

where * is the vectorial Hodge operator of the normal space, with orientation chosen so that r = — meaning 
that r is future-pointing for outward-pointing ^. Then r is normal to the horizon, g(^, r) = 0, _Lt — 0, has equal and 
opposite normalization g{T,T) = —g{£,,S,) and is regular in the null limit, becoming null itself, r — > ^. 

In spherical symmetry, the Kodama vector x provides a preferred flow of time, reducing to the stationary Killing 
vector for Schwarzschild and Reissner-Nordstrom black holes. It has a dual relation to the energy E which can be 
written as L^^E = Ar(x,T), for any normal vector ^ and its orthogonal dual r. This vector can be generalized by 

X = *idR)i (64) 

or the curl of R in the normal space, with components x — e^ {L^Rl_ — L^Rl^). Then x is orthogonal to R and 
the transverse surfaces, x ' dR = 0, _Lx = 0, has normalization ^(XjX) = —g^^{dR,dR) and becomes null on a 
trapping horizon, g{x,x) — Oi with x — i{dRy for 9± ^ 0. Flow lines of x and (di?)' are sketched in Fig. |21for 
typical gravitational collapse to a black hole, assuming cosmic censorship; for a comprehensively analyzed case, see 
Christodoulou 55] for the massless Klein-Gordon fleld in spherical symmetry. For a future outer trapping horizon 
[T^ . the area-radius vector (di?)" is spatial and x is temporal just outside the horizon, and (di?)' is temporal and x 
spatial just inside the horizon. 

In terms of these vectors, there is a remarkably simple and manifestly invariant expression for the matter energy 
density: 

e,„L^i? = r(x,T). (65) 

This holds for any foliation of spatial surfaces in any space-time, generated by a normal vector ^ with orthogonal dual 
T, without any gauge conditions. Then the integrated matter flux is 

[E],n = / ^Tix,T)Adx. (66) 
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FIG. 2: (colour online). Gravitational collapse: Penrose diagram of typical black-hole formation satisfying cosmic censorship, 
indicating the trapping horizon (blue bold line) , null infinity (straight lines) and the centre R — (curved line) , which is regular 
outside the trapping horizon and singular and spatial inside. Advanced time x'^ runs diagonally up-rightwards and retarded 
time x~ runs diagonally up-leftwards. Flow lines of the area-radius vector (di?)' (green short dashes) and the Killing-like vector 
X (magenta long dashes) are also indicated. Note the causal switch-over between x a-nd {dR)^ at the trapping horizon. 



Here the rescaling freedom x ^ x{x) in ^ and therefore r is cancelled by dx to leave an invariant expression. With 
the Ashtekar-Krishnan gauge conditions / = 0, = —1 and L^R = 1, one finds that x reduces to l+y^2/g^ on 

a trapping horizon, revealing that their permissible vector fields would generally coincide with x/L^R. The current 
formulation in terms of x is more manifestly invariant and physically interpretable, with x plf^ying the role of a 
stationary Killing vector. 

The integrated flux due to gravitational radiation can similarly be written as 

[E]g^ J *eix,T)Adx (67) 

where the effective gravitational-radiation energy tensor Q is defined to have components which are determined as 

e±± = lk±||V32^ (68) 
0±T = IC(±)lV8^e-^. (69) 

It has not been shown explicitly here that is a tensor in the normal space, since it is more easily seen in a spinorial 

formulation [s^]. Note that O is generally not symmetric; actually 8^ ^ Q |- with the admissible gauge choice 

/ = 0, but the general expressions will be retained here, in view of future generalizations. Since its components are 
non-negative, O satisfies the dominant energy condition, implying that the gravitational radiation carries positive 
energy. The four components of Q are interpreted as the energy densities of gravitational radiation: 0++ for ingoing 

transverse radiation, Q |_ for ingoing longitudinal radiation, for outgoing longitudinal radiation, and Q for 

outgoing transverse radiation, for = 0. Note also that (cr±, C(±)) each have the correct number (two) of independent 
components for describing the respective radiation. The identification and neat division of these four modes is another 
success for the dual-null method. 

The generator- volume form of the first law becomes 



[E]^ I ^{T{x,T) + eix,r))Adx. (70) 

In differential form, 

L^E^ (b*{T{x,T) + Q{x,T)). (71) 
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If the trapping horizon is spatial, one can use the unit normal vector f = t j ^g^x, recalling that gxx = 3(^1 = 
—g{T,T), to give the proper- volume form 

[E]= J i(T(x,f) + e(x,r)) (72) 

which is closest to that of Ashtekar & Krishnan, having identified O and x- These three forms perhaps most clearly 
demonstrate the nature of the first law as an energy-balance equation, expressing the increase in the mass-energy E 
of the black hole due to the energy densities of the infalling matter and gravitational radiation. 



IX. SURFACE GRAVITY AND A GIBBS-LIKE EQUATION 



In spherical symmetry, there is a natural definition of surface gravity k for a trapping horizon, satisfying an 
equation with the same form as the usual surface gravity for stationary black holes, but with the stationary Killing 
vector replaced by the Kodama vector p2, 56] . The energy-supply term / • "0 iii the first law can then be rewritten 
as kL(^A/8tt, yielding an equation with the same form as the first law of black-hole statics, which really is analogous 
to the Gibbs equation rather than the first law of thermodynamics. In seeking a general definition of surface gravity 
for trapping horizons, it is perhaps useful to note first that a formal first law can be given for any energy density w, 
if the surface gravity is related to it by 



An [ E - R (b *w] . (73) 



This yields the identity 



dE = ^-jl *KdA+ I *wdR + Rd(^] . (74) 



Since L^{E/R) = on a trapping horizon, projecting the above identity along the trapping horizon yields 

L^E ^ j) *kL^A+ (j) *wL^R. (75) 

This will be the desired Gibbs-like equation for black holes, generalizing that found in spherical symmetry plj. with 
a similar form to the first law of black-hole statics. 

One wishes to define surface gravity K(^±^ on a horizon with 6± = and work density W(^±^ as above. This determines 
K(^±) up to total divergences, and an argument in the next section fixes it as 

Ki±j^-jeU2L^0± + e+9.). (76) 

Here it is merely checked that K(^±) and u'(±) are related as above. The integrand in H73II can be rewritten using the 
cross-focusing equation (^5)) as 

e^ {2L^9± + = ef {\6ttT+^ - 0+0^) - 3? + 2|C(±)|^ - 2D ■ Cf±). (77) 

The last term integrates to zero, while the second and third terms constitute the integrand of the Hawking energy 
(EH), integrating to -IGnE/R. Then 

(j) = 4nE - AttR j * ^e^ r+_ + ^^^^^ = An (^E - R <j> *W(±)^ (78) 

as claimed. A version of the Gibbs-like equation H75|l was given previously in terms of the averaged surface gravity 
/*(«(+) + K(_))/2A and the averaged work density / *(iy(4_) + W(_))/2A |H3|- The first law resolves such 
ambiguities by determining the appropriate w and therefore k. The definition of outer trapping horizon |l4j . L^9± < 
(for 0± =0), also indicates the suitability of K(±) as measures of surface gravity, since K(-i-) > for such horizons, 
with K{±) vanishing somewhere on degenerate horizons. 

The standard definition of surface gravity for stationary black holes is via the formula (V Ax'') = kx' on a Killing 
horizon, where x is the stationary Killing vector and a flat (b) denotes the covariant dual (index lowering). The same 
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formula holds for the preferred time vector % on a trapping horizon, for an average surface gravity R, defined below. 
One finds = L_Rdx~ - L+Rdx'^ and 2V A = dx^ = 2L_L+i?dx+ A dx^ , noting that L+L_R = L^L+R since 
R is an integral scalar of the transverse surfaces. Then 

X • (V A x^) = RdR (79) 

where 

K = -e^L^L+R. (80) 
In particular, on a trapping horizon with 9± = 0, x' = ±dR and so 

X • (V A x') = ±nx^ (81) 

as desired. Thus R is guaranteed to recover the standard expression for stationary surface gravity if x reduces to the 
stationary Killing vector. If the expansions are constant on the transverse surfaces, Dd± = 0, then d± = 2L±R/R, 
so that K and K(-i-) all coincide. In particular, in spherical symmetry, k = R recovers n = 1/Am for Schwarzschild and 

K = ^J'rr^^~~(p- 1 (m + \J m? — for Reissner-Nordstrom black holes, where m is the mass and q the charge. The 
last result is non-trivial, and has not been recovered by other definitions of surface gravity for dynamical black holes. 

For Kerr black holes, a dual-null foliation giving the correct x and n is not known, despite the recent construction 
of a pair of dual-null foliations generating the horizons [s^. For the record, the results are as follows, using tildes 
to denote quantities associated with those particular dual-null foliations and reserving x = d/dt + 0.8/ d(f) and k, — 
\J m? — /2mr+ for the correct values, where m is the mass, ma the angular momentum, (i, r, 9, (p) are Boyer-Lindquist 
coordinates, r+ = m + \/ ni^ — a^ is the coordinate radius of the outer horizon and = a/2mr'+ is its angular velocity. 
The results are x — Qx ^nd k — Qk where Q = R?{dRI dr) /Tj, E = + cos^ and R = {r^ + a'^r^ + 2ma^r)^/'*. 
Thus, while x a-nd R, are out by a factor, it is the same factor. As a practical procedure, one could rescale x to x 
to obtain the correct k. It makes no difference to use R, as 9± = 2L±R/R in this case. This is clearly less than 
satisfactory, though there seems to be no other notion of surface gravity for dynamical black holes without some 
ambiguity or apparently ad-hoc procedure to recover the Kerr surface gravity. A guide here could be the issue of 
finding a dynamical zeroth law, stating that the surface gravity becomes constant as a trapping horizon becomes null. 
For R, one obtains DR ~ if Df = 0, which is a legitimate gauge choice on a trapping horizon of any signature. 
For the condition for = is D{e-^ L-9+) = 0, which constrains the null normal L generating a dual-null 

foliation away from a null trapping horizon. These issues are currently unresolved. 



X. A CLAUSIUS-LIKE EQUATION, ENTROPY FLUX AND ENTROPY CONSERVATION 

Given that stationary black holes theoretically possess an entropy A/ A, the parallel between the Gibbs-like equation 
(|75|l and the thermodynamic Gibbs equation ^ suggest defining a geometric entropy 

S ^ A/4 (82) 

for any trapping horizon |59l |. Units are such that the Planck and Boltzmann constants arc unity. The final issue to 
be addressed here is the analogue of the second law of thermodynamics (O , which for a non-isolated system states 
that the entropy of the system, minus the entropy supplied to the system, is non-decreasing; that is, entropy is either 
produced or conserved, but not destroyed. Then one needs a definition of geometric entropy supply. 

By comparing the first law (|5(J|I and the Gibbs-like equation (|75|l . it was shown indirectly that ^ • i/; = 
f *KL^A/8TrA. This is now shown explicitly and the result localized. Combining (|49I27I35I76|I successively yields 

^ ^ ^ P-^L,R ^ - '-^L,R ^ ^L,R (83) 



and so 



^e-V' = ^. (84) 
87r 



Thus the energy flux through the horizon is unexpectedly proportional to the surface gravity. Now recall as in §11 that 
heat flux q is classically proportional to temperature thereby defining entropy flux ip — g/i?, and that stationary 
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black holes possess a Hawking temperature "d = n/in. Then the above result suggests defining a geometric entropy 



lux 



(85) 

This argument is reminiscent of the original definition of entropy due to Clausius, who argued that the heat supplied 
to a system, divided by temperature, should be a total differential for closed cycles. Here is not the classical 
thermodynamic entropy flux of the matter, but should be regarded as a speculative definition of geometric entropy 
flux induced by matter (or gravitational radiation) in a strong gravitational field. In other words, if in some theory of 
quantum gravity one could indeed establish that black holes have a geometric entropy A/ A, perhaps one would also 
find that infalling energy fluxes induce corresponding geometric entropy fluxes. 
The corresponding geometric entropy supply So lO is given by 

L^So = <j>*S,-ip'^ L^A/A. (86) 

The second law of thermodynamics (jSJ for trapping horizons would state that L(^S > L^So- However, we have 
L^S ~ L^A/A and so 

L^S = LcSo. (87) 

Thus geometric entropy is conserved. This may sound radical, since what is normally called the second law for black 
holes is an inequality, like the second law of thermodynamics. A more faithful comparison with thermodynamics 
has shown that, while the geometric entropy of a black hole generally increases, it does so by the geometric entropy 
supplied to it by the infalling matter and gravitational radiation, with no net entropy production. This might perhaps 
be expected, since General Relativity is a classical theory which is symmetric under time reversal. Entropy can 
be produced by quantum-mechanical or statistical effects, in either the matter or quantum gravity, but appears to 
be absent classically. Similarly, the generalized second law should be stated as expressing combined (matter plus 
gravitational) entropy production, rather than entropy increase. 



XI. CONCLUSION 



The main results are summarized as follows, (i) An effective energy tensor O Ht)8Ht)9|l for gravitational radiation has 
been identified for dynamical black holes, (ii) A Killing- like vector x l|64f) . providing a preferred flow of time outside 
a dynamical black hole, has been identified and used to characterize the Ashtekar-Krishnan permissible vector fields, 
(iii) The Ashtekar-Krishnan energy-balance equation for dynamical black holes has been re-derived in a dual-null 
formalism, emphasizing geometrical invariance and corroborating the physical interpretation as an energy-balance 
law (|7UII72ll . with the black- hole mass-energy E growing due to the energy densities of the infalling matter and 
gravitational radiation, T(x,t) and 0(x,r). (iv) While the original proper-volume form of the law applies only to 
spatial (dynamical) horizons, here generator-volume and differential forms have been derived, which are both regular 
in the physically important limit of null (isolated) horizons, where a black hole is starved and ceases to grow. The 
new forms apply to any trapping horizon, thereby describing inner black-hole horizons, white holes, cosmological 
horizons, traversable wormhole mouths and evaporating black holes, (v) The energy terms have been divided into 
those which vanish if and only if the horizon is null, and those which generally do not (|50|l . and interpreted respectively 
as energy-supply and work terms, in analogy with the first law of thermodynamics (^3). (vi) A new definition of surface 
gravity k (|76f) has been given for dynamical black holes, such that the energy supply can be written in terms of k 
and area A, (|84() . just as in the so-called first law for stationary black holes, which is instead analogous to the Gibbs 
equation © of thermodynamics. An average surface gravity (|81|l has also been defined with respect to x by the same 
formula as stationary surface gravity, (vii) Since the energy-flux covector ip is proportional to k on a trapping horizon, 
the original Clausius concept of entropy suggests deflning an entropy flux 2T:ijj / k ()85f) , and it follows that entropy is 
conserved, ()87|l . for dynamical black holes. 

The last, perhaps surprising result does not contradict the fact that black holes grow, either by Hawking's area 
theorem for event horizons ^ ^ J,, ^^J or the area law for future outer trapping horizons , which instead reflect 
the fact that, classically, a black hole is the ultimate absorber. The black-hole area and presumed entropy A/A 
increase, but only by the entropy supplied to it by the infalling matter and gravitational radiation. It should be 
acknowledged that this physical interpretation is speculative, since it not known that dynamical black holes have a 
truly thermodynamical entropy A/ A and a local temperature K/27r, in the same way as is known for stationary black 
holes. It nevertheless illustrates that the classical first and second laws of black-hole mechanics la, Ifl Si are both 
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misnomers. While versions of both laws now exist for dynamical horizons, in addition, there are now more correct 
analogues for black holes of the first and second laws of thermodynamics: conservation laws for energy and entropy. 

The most practical result here is probably (i), in the context of current efforts to predict gravitational waveforms 
produced by dynamical black holes. It has generally been believed that gravitational radiation is well defined only 
in weak-field regimes, rather than the strong-field regime calculated by numerical simulations, thereby posing the 
radiation-extraction problem: how to extract the outgoing gravitational radiation and observable space-time strains 
from the numerical simulations. Now, however, the energy tensor defines gravitational radiation in the strong- 
field regime. Miraculously, a dynamical black hole itself provides the required structure, the dual-null foliation of 
ingoing and outgoing wavefronts generated from the trapping horizon, as located by existing numerical methods. 
The gravitational radiation may therefore be extracted from existing simulations by numerically implementing the 
coordinate transformation to the preferred dual-null foliation and calculating relevant quantities such as the conformal 
shear a-/R, which yields the Bondi news [13,1^, and the conformal energy flux R'^ip-. The variables are simply 
related to the actual strain tensor e/R to be measured by a gravitational- wave detector at large distance R, via the 
conformal strain tensor e — J{a-/2R)dx~ j23|. It should be noted that 6 is uniquely defined for a given trapping 
horizon. If one wishes to have a definition which is unique for a given space-time, an appropriate concept appears 
to be the trapping boundary [l4j . the boundary of an inextendible region whose every point lies on some trapped 
surface. Under smoothness assumptions, a trapping boundary is also a trapping horizon. However, in practice, 
trapping horizons have exactly the right level of uniqueness, since marginal surfaces are typically found numerically in 
any reasonable slicing of a black-hole space-time. In an asymptotically flat space-time, the waveforms must converge 
in advanced time and be equivalent for different slicings, so the practical issue is whether the waveforms converge 
sufhciently within the numerical domain; if so, it seems likely that they will converge to equivalent waveforms for any 
reasonable slicing. 

Strictly speaking, the results here have demonstrated the contribution of only the {0±±, 6±=p) components of the 
energy, with the (0::pzp, 0:p±) components determined by a symmetric treatment of trapping horizons with 0± = 0. 
However, the results generalize to uniformly expanding flows of the Hawking energy away from a black hole; the basic 
calculations may be found in |50j | , with the formally identical energy-tensor form to be reported subsequently '62] . 
An alternative (and probably easier) numerical implementation would therefore be to locate the level surfaces of 
these flows in the original spatial hypersurfaces, calculating desired quantities separately for each hypersurface. The 
Hawking energy can also be used to recover the Bondi energy loss at future null infinity [sol l5lj , which can also be 
written in a formally identical energy-tensor form. Thus Q apparently provides a good physical measure of the energy 
densities of gravitational radiation all the way from a black-hole horizon out to infinity. The radiation-extraction 
problem for dynamical black holes has thereby been theoretically solved. 
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